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Two-particle two-hole contributions to electromagnetic response functions are computed in a 
fully relativistic Fermi gas model. All one-pion exchange diagrams that contribute to the scattering 
amplitude in perturbation theory are considered, including terms for pionic correlations and meson- 
exchange currents (MEC). The pionic correlation terms diverge in an infinite system and thus are 
regularized by modification of the nucleon propagator in the medium to take into account the finite 
size of the nucleus. The pionic correlation contributions are found to be of the same order of 
magnitude as the MEC. 
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I. INTRODUCTION 

The goal of this paper is to present a fully relativistic 
calculation of the two-particle two-hole (2p-2h) contri- 
butions to the inclusive (e, e') response functions of nu- 
clei for intermediate to high momentum transfers in a 
Fermi gas model. Consistency with perturbation theory 
is maintained and all diagrams with one-pion exchange 
in the nuclear current are considered, constructed by at- 
taching a photon to all possible lines in the basic one-pion 
exchange Fcynman diagram. In this way not only meson- 
exchange currents (MEC) arise (for example, where the 
photon is attached to the pion), but also pionic correla- 
tion diagrams, where the virtual photon is absorbed by 
one of the two interacting nucleons. Both kinds of di- 
agrams are considered in our model, together with the 
usual virtual A-isobar electroexcitation and decay. 

We are motivated by previous work presented in [HQ, 
where only the MEC were included in the 2p-2h trans- 
verse (T) response, together with earlier work both in 
non-relativistic Q and relativistic 0-0] regimes. The 
contribution found from the 2p-2h excitations is small 
at the quasielastic (QE) peak, and increases with energy 
transfer, being more important in the dip region, where 
it is dominated by the A current. At the non-relativistic 
level attempts were also made to evaluate the 2p-2h con- 
tribution of MEC in the T-response for finite nuclei in a 
shell model 

The MEC are not the only two-body operators able 
to induce 2p-2h excitations. The correlation operators 
arising from Feynman diagrams where the photon is at- 
tached to a nucleon line, exchanging a pion with another 
nucleon, are of the same order as the MEC in the per- 
turbative expansion and should be included to be con- 
sistent |10l4l2| . These diagrams, however, present the 
problem of giving an infinite answer in a Fermi gas model. 



The reason is that there is a nucleon propagator that can 
be on-shell in the region of the quasielastic peak. Since 
the response function is the square of the amplitude, the 
resulting double pole gives an infinite result after integra- 
tion. In dealing with this problem, in [T^] a prescription 
was followed by keeping the lines with a nucleon propa- 
gator strictly off the mass shell. A different approach was 
taken in by subtracting from the proper self-energy 
its value on the mass shell, with the unphysical shortcom- 
ing of obtaining negative results for the 2 p-2h responses 
to the left of the QE peak. Finally, in [12( a nucleon 
self-energy in the medium was introduced in the nucleon 
propagator. In dealing with the seven-dimensional inte- 
grals appearing in the 2p-2h responses, some of the pre- 
vious calculations have resorted to the approximation of 
setting the two hole momenta both equal to zero in some 
of the diagrams flfj| or by taking into account only an 
average nucleon momentum (l2| . 

In this work we revisit the double-pole problem to an- 
alyze the nature of the divergence of the resulting contri- 
butions. By isolating the divergent terms we are able to 
link them to the infinite extension of the Fermi gas sys- 
tem. In fact the double-pole term can be related to the 
probability of one-nucleon emission followed by nucleon 
re-scattering off another nucleon, with the final ejection 
of two particles. This probability is infinite, since it is 
proportional to the propagation time of a real nucleon 
in a Fermi gas. This fact was pointed out in [l2| where 
it was cured, as mentioned above, by introducing a nu- 
cleon self-energy with an imaginary part giving it a finite 
lifetime for collisions. In this paper we use a similar pro- 
cedure by introducing a finite imaginary part it in the 
nucleon propagator, but with a new meaning for the free 
parameter e. Instead of being an imaginary part of the 
nucleon self-energy for collisions, we relate it to the time 
T that a nucleon can travel across the nucleus before 
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leaving it. Hence this term accounts for the finite size of 
a real nucleus in contrast to an infinite system like the 
Fermi gas, where T is infinite. The value of e can be esti- 
mated to be roughly about 200 MeV, appreciably larger 
than the usual values of the nucleon width for collisions. 

The structure of this work is as follows. In Sect. II we 
present our model and define the 2p-2h response func- 
tions and the two-body current operators. We discuss in 
depth the divergence of the correlation diagrams and the 
need to introduce the parameter e in Sect. Ill (details of 
the numerical calculation are given in the appendices). 
In Sect. IV we present results for the 2p-2h longitudi- 
nal and transverse response functions. In the case of the 
correlation diagrams we present results for several values 
of the parameter e. Finally, in Sect. V we present our 
conclusions. 



II. MODEL FOR 2P-2H RESPONSE 
FUNCTIONS 

We consider an electron that scatters off a nucleus 
transferring four- momentum = (w,q), with u the 
energy transfer and q the momentum transfer. We fol- 
low closely the notation of [13| . Assuming plane waves 
for the electron, working in the laboratory system and 
taking the z direction along the momentum transfer, the 
inclusive cross section is written as 



da 



dQLdu) 



<*M [v L RL(q,u) + v T R T (q, to)} , (1) 



where cjm is the Mott cross section, vl and are the lep- 
ton kinematic factors, and the relevant quantities are the 
longitudinal RL(q,uj) and transverse Rp(q,Lo) response 
functions, respectively. These are defined as the follow- 
ing components of the hadronic tensor, 

W 00 
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(2) 
(3) 



where 



Y,{f\^{Q)W{f\J v iQ)W{Ei + u- E f ) (4) 



and J M (Q) is the nuclear current operator. 

In this paper we take the initial nuclear state as the rel- 
ativistic Fermi gas (RFG) model ground state, \i) = \F), 
with all states with momenta below the Fermi momentum 
kp occupied. The sum over final states can be decom- 
posed as the sum of one-particle one-hole (lp-lh) plus 
two-particle two-hole (2p-2h) excitations plus additional 
channels. In the impulse approximation the lp-lh chan- 
nel gives the well-known response functions of the RFG. 
Here we focus on the 2p-2h channel where the final states 
are of the type |/) = p'lS^, p' 2 s' 2 , h^si, h 2 " 1 s 2 ), where 
are momenta of relativistic final nucleons above the 
Fermi sea, p\ > kp, with four- momenta P[ = (£^,p£), 
and Hi = hj) are the four-momenta of the hole states 
with hi < kp- The spin indices are and Sj. 



A. 2p-2h Response functions 

Since we have two species of nucleons, the 2p-2h 
responses can be further decomposed as the sum of 
two-proton (PP), two- neutron (NN) and proton- neutron 
(PN) emission 



Rk = Rk(PP) + Rk(NN) + R K {PN) 



(5) 



For the PP channel we write down the L response as 
(likewise for the T response): 

Rl(PP) = 

JEEEE |(pip 2 hr 1 h 2 -v o (Q)i^r 

p'jSj p' 2 s' 2 hisi h 2 s 2 

x 5{E[ + E' 2 - ui - E l - E 2 ) , (6) 

where the spin indices are implicit in the matrix elements. 
The factor j comes from anti-symmetry of the wave func- 
tions, to avoid double counting of the final states under 
the interchange of the indices 1' <-> 2' and 1 o 2. Ex- 
ploiting the anti-symmetry, the many-body matrix ele- 
ment of a two-body operator can be written as the direct 
minus exchange part of the two-body current matrix el- 
ement 

(pip^hrV = (p , 1 p 2 |J M |h 1 h 2 )-(pip 2 |j^|h 2 h 1 ) 

which we write in terms of the two-body current function 
j M (p'i! P2j hii n 2) to be specified below, 

(p' lP2 1 J"|hih 2 ) = (2 7 r) 3 <5(p' 1 + p 2 - hi - h 2 - q) 



Tfl 



(7) 



Going to the thermodynamic limit and integrating over 
the momentum p' 2 wc obtain 



Rl(PP) 



V 



E 



dVi d a hi d A h 2 
(2tt) 3 (2tt) 3 (2tt) 3 



'FTnv |i (Pl>P2, h l, h 2) 



EiE2E' 1 E' 2 
x 5(E[ +E' 2 -u>- El - E 2 )6(p 2 - k F ) , 



(8) 



where p 2 = hi + h 2 + q p^ , and the integration limits 
are hi, hi < kp, p[ > kp. We have defined the anti- 
symmetrized current function 

f(Y, 2', 1, 2) A = f{V, 2', 1, 2) - 2', 2, 1) 

with obvious meaning for the abbreviated arguments. 
Expanding the square inside the integral in Eq. © , three 
terms are obtained: 

|^(1',2',1,2) A | 2 = \f{l',2',\,2)\ 2 + \f{l',2',2,l)\ 2 



2Re^(l',2',2,l)*^(l',2',l,2). 



(9) 
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FIG. 1: MEC diagrams considered in the present study. Di- 
agrams (a,b) correspond to the seagull, (c) to the pionic, and 
(d-g) to the A current, respectively. 



Changing variables 1 f> 2 in the second term under the 
integral, we obtain the first term again. Hence we can 
finally write for the PP response 



/? (p P ) = I V f rfVl d " kl d3h2 
M ' ~ 2 ^ J (2tt) 3 (2tt) 3 (2tt) 3 



EiE 2 E' l E' 2 
-Re j°(p'i>P2, hi, h 2 )*i (p'i,P2, h 2 , h x )] 



x S(E[ + E 2 



E 1 -E 2 )9(p' 2 -k F ). (10) 



Note that the factor \ in front of the sum comes from the 
anti-symmetry of the particles (protons). A similar ex- 
pression is obtained for the NN response Rl(NN). In the 
case of the PN channel we subtract the charge exchange 
contribution without any symmetry term because there 
are no additional isospin sums, and the result is 



R L (PN) = V ]T J 



d 3 p[ d 3 /ii d 3 h 2 
(2^ (2^)3 (2^)3 



rrr 



— KPNlfip'^PzMM^PN) 



Ei E2 E^ E'2 
x 5(E[ +E' 2 -uj-Ei- E 2 )6(p' 2 - k F ) 



[11) 



Finally, note that the 2p-2h response is proportional 
to the volume of the system V which is related to the 
number of particles M (protons or neutrons) by V = 
3n 2 Af/k 3 F . 



B. Two-body current matrix elements 

The MEC considered in this work are represented 
by the Feynman diagrams of Fig. 1. The pionic four- 
momenta K\ , K 2 are defined via K'f = P'^ — as the 



FIG. 2: Correlation diagrams considered in the present study. 
Diagrams (a,b) correspond to the forward, and (c-d) backward 
contributions, respectively. 



four- momenta given to the nucleons 1 and 2, respectively, 
by the exchanged pion. 

Assuming pseudo-vector nuclcon-pion coupling, the 
fully relativistic two-body current matrix elements are 
given by [H[l|: 



(a-b) Seagull or contact: 



^(Pl,P2iPl>P2) = ^r«e3a6"(Pi)T- a 75 ^iw(pi) 



FY 



K\ — ml 



u(p' 2 )nj 5 ^u{p 2 ) + {1 ^ 2) 



112) 



(c) Pion-in-flight: 



ip(Pi,P2>Pl>P2) 



F w (Kx - K 2 y 



ml le3ab (K?-ml)(Kl-ml) 



XM(pi)r a 7 5 ^Tiu(pi)u(p 2 )r fc 7 5 fi 2 u{p 2 ) 



(13) 



In the above we use the Einstein convention for the sum 
over a repeated isospin index a = 1,2,3. Moreover, F^ 
and F n are the electromagnetic isovector nucleon and 
pion form factors, respectively. The spinors are normal- 
ized according to the Bjorken and Drell convention [l5| 
and the pion-nucleon coupling constant is f" /Att = 0.08. 

• (d-g) Delta current: 



j£(Pi,P 2 ,Pl,P2) = ^It 1 ^ 2 ^' l )T^l)u(p 1 ) 



m\ K 2 — m\ 



Xu(p 2 )-r 75 ^2M(P2) + (1 o 2) . 



(14) 



The vector T£(l) is related to the pion electroproduction 
amplitude 

W) = K 2ja Q a ^G^ p {H 1 + Q)Sf{H l )T a Tl 

+T z TiS£ p {P[)G%{P[ - Q)e*»tf 2 , a .(15) 

The forward A electroexcitation tensor is [24| 

Sf{H x ) = G" 4 [gi Q - g 2 H x ■ Q + g 3 Q 2 ] 75 

- e^Q v [ gi ^-g 2 H^+g 3 Q^ 5 (16) 
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and the backward tensor amplitude is 

sr(P{) = 75 [gi Q-nPi-Q- 93Q 2 ] &*" 

- 75 [gir - 92P'S - 9zQ»\ Q v Q vp ■ (17) 



The tensor M „ is defined by 



= 9^u - -.lulu ■ 



(18) 



For the A propagator we use the usual Rarita-Schwinger 
(RS) tensor 



GfJP) 



P + m A 

p2 



'A 



X 



gp? -pel? -3-^ 



2 PpP p lp p p - 7p P f 



(19) 



o A 3mA 

In what follows we perform the substitution to a — > toa + 
|r(P) in the denominator of the propagator to account 
for the A decay probability. Finally, the electromagnetic 
coupling constants gi are given by 



gi 



2niN 



gi 



G 2 
Am 2 N 



.93 



G 3 
Am 2 N 



(20) 



Our approach for the A follows, as a particular case, 
from the more general form of the -fNA Lagrangian of 
Pascalutsa et al. [l6[. The A coupling constants used 
here are Gi = 4.2, G 2 = 4, G 3 = 1, and / iJV a = 4x0.564. 
The correlation current is defined in Fig. 2, and given 

by 



Jcor(PiiP2>Pl)P2) 



p _ 1 

— u(p' 1 )r a75 ^iu( Pl )— 5 2 

x u(p 2 ) [r a75 ^i5 F (P 2 + Q)r"(Q) 

+r^(Q)S F (^ - Q)r a75 fo] u(p 3 ) 
+ (lo2), (21) 

where Sf (P) is the Feynman propagator for the nucleon 



S F (P) 



p2 _ m 2 _|_ j £ 

and T' i ((5) is the electromagnetic nucleon vertex, 



r**(Q) = Fx-f + — F 2 ^ V Q U 
2m 



(22) 



(23) 



The nucleon form factors F\ and F 2 are given by the 
Galster parametrization [T^j- 

The isospin sums and isospin matrix elements must be 
performed separately for each isospin channel. Explicit 
expressions are given in Appendix A. 



III. DIVERGENCE OF THE CORRELATION 
RESPONSES 

The response functions computed using the correlation 
current in Eq. (|21[) arc divergent in the Fermi gas. There 



are two sources for this divergence: the first one comes 
from the double pole of the propagator when taking the 
square of the current. This divergence can be shown to 
behave as 1/e plus principal value terms going as loge. 
The second source is related to the behavior of the prin- 
cipal values arising from the double and single poles near 
the RFG boundary of the quasielastic peak, where the 
principal values present a logarithmic divergence. 

To illustrate the mathematical structure of this diver- 
gence we isolate as an example the singularities produced 
by the diagram of Fig. 2(a). The corresponding current 
operator can be written as 



3 



Ei 



£1 



hi+q 



(24) 



where E p = m? + p 2 is the on-shell energy. We have 
explicitly extracted the divergent part of the denomina- 
tor, with a pole for 



Eht+q — El 



(25) 



in the limit e — > 0. The above equation is equivalent 
to the quasielastic condition for emission of an on-shell 
nucleon with four-momentum Hi + Q. In fact, for a given 
value of hi, Eq. (|25[) holds when the angle between hi 
and q is given by 



cos 61 



Q 2 + 2Eiuj 
2h x q 



(26) 



Since the condition — 1 < cos#i < 1 defines the boundary 
of the quasielastic peak, the pole can always be reached 
in that region. 

To study the behavior of the response functions due to 
this pole, it is convenient to change the variable 81 to a 
new variable defined by 



Xi — Ei + uj — E] 



hi+q 



(27) 

in the integral over hi in Eq. (|10|) . Then the components 
of the total current matrix element can be written as a 
function of xi in the general form 



Xi + ie 



+ g(xi) , 



(28) 



where the first term comes from diagram 2(a) and the 
function g{xi) comes from the sum of the remaining dia- 
grams, and is finite for xi = 0. Since the current appears 
squared in the response function, we are dealing with the 
integral of a function of the kind 



— -2+\g(xi)\ +2Re — — — . (29) 



x 



xi — ie 



When integrating this function over x\, and taking the 
limit e — > 0, the first term has a double pole for xi = 0, 
while the third one has a single pole. To deal with the 
single pole we use the usual Plemcli relation, 



1 



x + ie 



V iir5(x) . 

x 



(30) 



5 



To apply a similar relation for the double pole term, we 
add and subtract the on-shell value \ip{Q)\ 2 /{x\ + e 2 ). 
Taking the limit e — > we can use relations which are 
valid for any function ip( x ) 



4>(x) - <£(0) 



dx 



^)-V(O) 



dx (31) 



and 



Then Eq. 



dx 



tan — h tan 
e 



i>(0) ~ -v(o). 

(32) 



can be written in the form 



l/(*i)| 2 = V 



\<f(xi)\ 



m 



Re (jJ*(xi )q(xi ) „ . . , , „, 

2p g v lm i; _ 27rIm v »(o) 5 (o)J(xi 



-ttS(xi 



(33) 



The last 0(1/ e) term in Eq. (|33p provides the dominant 
contribution to the response function, being infinite for 
e — >• 0. Due to the <5 function, that term does not con- 
tribute outside the quasielastic-peak region, where x\ is 
different from zero. 

The principal values present in Eq. ([3"3"| also diverge in 
the particular case in which one of the limits of integra- 
tion is zero. In that case the principal value in Eq. (|30p 
should be computed instead using 



V 



ip(x) 



dx = 



i/>(x) - V(0) 



1 



dx + -ip(0) In ■ 



+ e' 
(34) 

and it gives a In e term if a or b is zero. That situation 
in fact occurs throughout the quasiclastic region, and 
in particular at the boundary of the quasielastic peak. 
Therefore one expects an additional divergence ~ 0(ln e). 

The meaning of the term nS(xi) is explained in 

what follows. Diagram 2(a), when the intermediate nu- 
cleon is on shell, gives the probability of a lp-lh elec- 
troexcitation times the probability of quasiclastic nucleon 
scattering. Since the interaction probability is propor- 
tional to the interaction time T, the probability of this 
re-scattering process is proportional to T 2 . Therefore the 
cross section is proportional to T. In an infinite system 
such as the Fermi gas, the intermediate nucleon never 
leaves the nucleus and therefore T — > oo. However, in a 
finite nucleus one expects no divergence because a high- 
energy nucleon will leave the nucleus in a finite time. 
Therefore the interaction time is finite. 

The relation between e and T can also be obtained by 
inspection of the momentum-space propagator in quan- 
tum field theory [l8|, computed as the vacuum expecta- 
tion value of time-ordered Fermion fields. The value e 
in the denominator of the propagator can be seen as a 
rcgularization parameter in the Fourier transform of the 



time step function for a particle with four-momentum 
P M = (po,p) 



T/2 

dte^-v^eit) = — 

-T/2 PO 



(35) 



where T — > oo and e — > 0. For a real particle, po—E p = 0, 
the left-hand side of the above equation is T/2, and the 
right-hand side is 1/e. Therefore 



T 
2" 



(36) 



This can be obtained alternatively by replacing the on- 
shell value of the propagator in Eq. (|35|) as a delta func- 
tion 



— = lim — 

e po->-E p p 



tt£(0) 



and using the integral representation 



5(0) 



1 



r T/2 rp 

lim — / dt = — . 



(37) 



(38) 



In this paper we cure the divergence of the correlation di- 
agram by a rcgularization procedure, using a finite value 
for e to account for the finite propagation time of a high- 
energy nucleon in a nucleus before leaving it. To estimate 
the value of e for a nucleus such as 12 C, we assume that 
the nucleon moves at the velocity of light and it has to 
cross a distance equal to the nuclear radius R ~ 2 fm. 
Then 



2h 



2hc 



400 



MeV~ 200 MeV. 



(39) 



Note that this value, e ~ 200 MeV, is very different from 
the nucleon width r ~ 10 MeV which is usually obtained 
in nuclear matter as the width for nuclear inelastic in- 
teraction. In practice the value of e can be taken as a 
parameter to be fitted to data. In the next section we 
perform a study of the dependence of our results upon e. 
Unless otherwise specified we assume e = 200 MeV. 

At this point we should mention that the use of 
Eq. (|29|) to compute the 2p-2h response functions be- 
comes impractical due to complications in the numerical 
calculation of principal values in multidimensional inte- 
grals including the four diagrams of Fig. 2 (and the cor- 
responding exchange parts). Since we are forced to use 
a finite value of e, it becomes more convenient to keep 
from the beginning the ie term in the denominator of the 
nucleon propagator in Eq. (|2"Tj) . 



IV. RESULTS 

Here we present results for the longitudinal and trans- 
verse response functions for inclusive two-particle emis- 
sion. We compute the 2p-2h response functions in 
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FIG. 3: 2p-2h transverse response of 56 Fe at q = 550 MeV/c. 
Three values of the parameter e are shown. Thin solid lines: 
Correlation only. Dotted lines: MEC only. Thick solid lines: 
total. Dashed: RFG OB results. 



the RFG model as the 9-dimensional integrals given by 
Eqs. (|10lll[) . The energy delta function can be used to 
integrate over p[ , fixing the energy E[ of the first parti- 
cle. More details are given in Appendix B. By rotational 
invariance considerations, one of the the azimuthal an- 
gles can be fixed, multiplying at the same time the result 
by a factor 2ir. We choose <f>[ = 0. At the end we have 
a 7-dimcnsional integration to be performed numerically. 
The usual procedure is to use a multi-dimensional Mon- 
tccarlo (MC) integration. Since the pole structure of the 
integrand is numerically delicate, in this work we use 
instead a mixed Montecarlo-Simpson integration proce- 
dure. The Simpson algorithm is used for integration over 
the angles of the two holes 6\ , 62 and of the first par- 
ticle Q' x . The remaining 4-dimensional integral over the 
hole momenta h\ , /12 and their angles <j>\ , <p2 is made by 
Montecarlo. To keep the CPU times manageable we use 
a number of MC points of the order of 10 3 for q = 1 
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FIG. 4: As for Fig. 3, but now at q = 1140 MeV/c. 



GeV/c. For other values of the momentum transfer the 
number of MC points is modified linearly with q. We have 
performed a study of the stability of the results with the 
number of MC points and have found that the error from 
the integration procedure is within a few percent. 

A pion-nuclcon form factor is included in the 2-body 
currents: F nNN (K v ) = (A 2 - m 2 )/(A 2 - A" 2 ), with A = 
1.3 GcV. We use the same value for the 7rAA form factor 
in the Delta current. The electromagnetic form factors 
are those of Galster for the nucleon, and those used in [l3|, 
[13 for the MEC. 

To make contact with previous work, wc apply our 
model to compute the 2p-2h longitudinal and transverse 
response functions for the nucleus 56 Fe, and for momenta 
q = 550 and 1140 MeV/c. The results are presented in 
Figs. 3-6, where the separate contributions of the corre- 
lation and MEC currents to the 2p-2h responses are also 
shown. The lp-lh responses produced by the one-body 
(OB) current in the relativistic Fermi gas (RFG) without 
interaction are also shown. 

A critical input for our model is the value of the pa- 




ramctcr e in the nucleon propagator, introduced to cure 
the divergence of the double pole. To see how the re- 
sponses for the correlation contribution depend on e we 
show results for three different values: e = 100, 200 and 
300 MeV. For e = 100 MeV, the correlation 2p-2h con- 
tribution presents a shape with a maximum in the region 
of the quasielastic peak, but with a long tail extended 
to high transferred energies. The maximum is reminis- 
cent of the pole structure of the nucleon propagator, and 
therefore a resonance appears for kinematics correspond- 
ing to the quasielastic condition in Eq. (f2"B|) . A shift to 
higher energies (of the order of ~ 40 MeV) is seen in the 
case of q = 550 MeV/c (Figs. 3, 5). Indeed for this value 
of q the phase space for two-particle emission causes a 
suppression of the low-energy side of the response func- 
tion. 

The resonant structure produced by the 2p-2h correla- 
tion contribution diminishes significantly with increasing 
values of the parameter e. Notice that for e > 200 MeV 
there is no maximum located at the QE peak. 

For an even lower value of the escape width, say e = 



50 MeV, the magnitude of the resonant peak is of the 
same size as the OB response function. This correction 
coming from 2p-2h states is obviously too large to be 
compatible with experimental data that are already of 
the order of the lp-lh response at the region of the QE 
peak. It should be mentioned that, although the 2p-2h 
contribution should be added to the lp-lh one, the latter 
should be first corrected for final-state interaction (FSI) 
contributions not included in the bare RFG results shown 
in the figures. In fact FSI contribute importantly to one- 
nucleon emission through the coupling of lp-lh to 2p-2h 
states in the final nucleus (liij ). These processes involve, 
in particular, two-pion exchange, and are therefore of 
the same order as the 2p-2h response in the perturbative 
series, since it is the square of one-pion exchange matrix 
element. The inclusion of such contributions is out of the 
scope of the present study. 

The dependence of the correlation responses on the pa- 
rameter e is better appreciated in Figs. 7 and 8, where 
we show its contribution for the three chosen values of 
e in the same plot. In the QE region the height of the 
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responses approximately reduces to one half when e dou- 
bles. This behavior follows because of the leading 1/e 
dependence in Eq. (|33|) . coming from the pole in the prop- 
agator. For high lu the results are more similar and they 
are almost independent of e in the high-energy tail. In 
this case, i.e. large lu, there is no pole in the integrand 
and the contribution from the propagator is less sensitive 
to the precise value of e. 
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FIG. 7: 2p-2h correlation contribution to the L and F re- 
sponses of 56 Fe for q = 550 MeV/c. Fhree values of the 
parameter e are shown. With dotted lines from up to down, 
e = 100, 200, 300, respectively. Solid lines: RFG one-body 
responses. 
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FIG. 8: As for Fig. 7, but now at q = 1140 MeV/c. 



Let us return now to Figs. 3-6, where the MEC sepa- 
rate contribution is also shown. The transverse response 
(Figs. 3, 4) has a large peak with a maximum around 
lu = (m 2 A + q 2 ) 1 ! 2 — 77ijv, that comes from the A propa- 
gator appearing in the A-current. It has the same reso- 
nant structure as the correlation current, but located in 
the region of the A peak, where the real pion emission 
cross section has a maximum. We do not include the 
pion emission channel in our calculation. Both channels 
should be summed up to obtain the total inclusive (e, e') 
cross section. 

The A peak is very small in the longitudinal response 
presented in Figs. 5 and 6. This is consistent with the 
predominant transverse character of the A current, hence 
providing a small contribution to the longitudinal chan- 
nel. For q = 550 MeV the MEC 2p-2h contribution 
is large (small) in the T (L) response. However, for 
q = 1140 McV (Figs. 4, 6) we find a larger effect in Rl 
coming from the MEC seagull and pionic at large energy 
transfer. Indeed in a non-relativistic expansion in powers 
of q/rriN the time component of the MEC is of higher or- 
der than the transverse one. However, for q = 1140 MeV, 
q/rriN is larger than one, and the relative L component 
of the MEC, compared to the T one, starts to increase. 

In the case of the correlation current, we observe that 
its contribution, compared with the OB responses, is sim- 
ilar in the T and L channels. Note that in the correlation 
current (Fig. 2) the photon couples directly to a nucleon 
with the same interaction vertex as the OB current. 
The other side of the diagram with a pion coupled to a 
second nucleon is independent of the particular compo- 
nent of the current. 

The separate effects of the different currents con- 
tributing to the 2p-2h transverse responses are shown in 
Figs. 9, 10. As shown, the seagull plus pionic (SPP) cur- 
rents alone give a small effect compared with the contri- 
butions from the A and correlations. In fact, for e = 200 
MeV the correlation response is much larger (by a factor 
2 or 3) than the SPP response function (middle panels 
in Figs. 9, 10). We also observe that the separate seagull 
contribution is larger in magnitude than the pionic one, 
which is negligible for q = 1140 MeV/c. Note that the 
two currents interfere destructively and partially cancel 
when both are considered in the SPP responses (bottom 
panels). 

In Fig. 11 we show the transverse response obtained 
by adding the total 2p-2h contributions to the OB re- 
sponse. A word of caution should be raised when analyz- 
ing these results. First, we have not added the correlation 
nor MEC corrections to the lp-lh response. Moreover, 
the two-pion-exchange interaction generates self-energy 
corrections to the OB current that lead to interference 
effects of the same order in the expansion as the correc- 
tions included here. As an example, FSI are known to 
redistribute the strength of the responses, producing a 
hardening, a reduction of the maximum and an increase 
of the high-energy tail [2(|. Recently also a large effect 
from both MEC and FSI has been found in the lp-lh 
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FIG. 9: Contributions to the transverse response of 56 Fe for 
q = 550 MeV/c. The dashed lines are the lp-lh response 
with OB current only. The rest of the lines are 2p-2h con- 
tributions, (a) Upper panel: Thin solid: Correlation only. 
Dotted: MEC only. Thick solid: total, (b) Middle panel: 
Thin solid: Correlation only. Dotted: seagull+pionic only. 
Thick solid: A only, (c) Bottom panel. Thin solid: pion-in- 
flight only. Dotted: seagull+pionic only. Thick solid: seagull 
only. 



channel for high momentum transfer 21] . which should 
be added to the present results. Finally, the process of 
real pion emission (not included here) gives also a con- 
tribution in the transverse response located mainly the 
region of the delta peak. 

So far we have presented results for intermediate to 
high momentum transfer. Results for lower values of 
q = 370 and 410 MeV/c are shown in Figs. 12 and 13 
for the T and L response functions. This allows us to 
compare the present results with previous non-relativistic 
calculations [l(J. In Fig. 12 the structure function 
St = t^-Rt is presented, to allow a direct compari- 
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FIG. 10: As for Fig. 9, but now at q = 1140 MeV/c. 



son with Figs. 11 and 12 of The separate MEC 

and correlation contributions to the 2p-2h T response 
shown in Fig. 12 arc similar to the ones presented in 
[To[ | . The MEC produces a tail above the QE peak that 
increases with the energy transfer. The presence of cor- 
relations lead to an additional, significant raise of the 
tail. Note that our correlation results are obtained for 
e = 200 MeV. In [l(| another prescription to deal with 
the nucleon pole was adopted. From our results we con- 
clude that both prescriptions are compatible numerically. 
The OB response of [l(| included RPA correlations pro- 
ducing a reduction and hardening of the OB response. 
The 2p-2h longitudinal responses were not computed in 
[Tol | , since the time components of the MEC arc of higher 
order in the non-relativistic reduction and hence, they 
were expected to be very small. However, our prediction 
for the correlation 2p-2h contribution in the L response, 
presented in Fig. 13, shows a similar effect as in the T 
response, i.e., & tail also appears for high energy trans- 
fer in the L response coming from correlations. Contrary 
to the T channel, MEC give no contribution in the L 




FIG. 11: Transverse response of 56 Fe at q — 550 and 1140 
MeV/c. Thin solid: Correlation only for e = 200 MeV. Dot- 
ted: MEC only. Thick solid: total one- plus two- body re- 
sponses. Dashed: RFG lp-lh response with OB current only. 




FIG. 12: Transverse structure function St of 56 Fe at q = 370 
and 410 MeV/c. The parameter e = 200 MeV. To compare 
with Figs. 11 and 12 of [l(J, St is defined as ^ S T = ^Rt- 
Thick lines: RFG lp-lh results. Dashed lines: 2p-2h, MEC 
only. Thin solid lines: 2p-2h total, MEC plus correlations. 
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FIG. 13: As for Fig. 12, but now for the longitudinal struc- 
ture function Sl = ^^Rl- 



response. 

Since the 2p-2h excitation is produced in this work by 
one-pion exchange, the results are strongly dependent on 
the details of this particular interaction. This is illus- 
trated in Fig. 14 where we show how the results depend 
on the strong irNN form factor for q = 1140 MeV/c. The 
results without a form factor, i.e., with F^nn = 1, are 
about three times as large as the results with the form 
factor. This is different from the findings at low momen- 
tum transfer Q , where the pion form factor can be safely 
ignored. 

Another issue is the dependence of the results on the A 
form factors used in this work, both the electromagnetic 
and the strong ones, which are somewhat different from 
theparametrization used in the 2p-2h MEC calculation 
of [l| . Calculations done with both sets of parameters are 
compared in Fig. 15. Our calculation gives a larger con- 
tribution for the T response than the one of pj. Hence 
the use of the same form factors reduces the discrepancy 
between the two calculations. Some of the remaining dif- 
ferences could be linked to other details of the models, in 
particular to the different Lagrangian chosen for A elec- 
troexcitation. We should note that the two models are 
fully independent. While all the spin sums are performed 
analytically in [l| resulting in thousands of terms to be 
numerically integrated, in this work we first compute the 
spin matrix elements of the current and later we evaluate 
the squares and perform the sums numerically. 

Before concluding, we would like to stress that the 2p- 
2h responses in the present model are crucially dependent 
on details of the pion interaction. A critical ingredient 
of the model is the value of the parameter e, identified 
with an escape width of a high-energy nucleon from the 
nucleus. We have proven that a value e ~ 200 MeV 
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FIG. 14: 2p-2h T response of 56 Fe at q = 1140 MeV/c. Thin 
lines with nNN form factor. Thick lines without form factor 
{F wNN = 1). (a)Total, (b) Correlations only, (c) MEC only. 



leads to results in agreement with the previous calcula- 
tion of |Toj]. This parameter e is different from the usual 
interaction width of particle states, usually associated 
with matrix elements of the phenomenological imagi nary 
optical potential derived from elastic scattering data [19| . 
It has also been computed in nuclear matter in a semi- 
phenomenological approach [23[ . The resulting width for 
100 MeV nucleons is of the order of 10 MeV, which is too 
small to give reasonable results in our calculation. This is 
due to the 1/e behavior of the 2p-2h response divergence 
in the QE region, where the pole is being hit. 

Due to this divergent behaviour, for e = 5 MeV the re- 
sults would be almost one order of magnitude larger than 
the OB responses at the maximum. We have checked that 
the the 1/e term in the forward diagrams is the main con- 
tribution to the 2p-2h correlations in the QE region for 
e > 20 MeV. 

The importance of correlations, for the same value of 
e, increases with the nuclear mass. We have checked that 
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FIG. 15: 2p2-h MEC only contribution to the transverse re- 
sponse of 56 Fe. Solid lines: computed with the A form factors 
used by [14| . Dashed lines: computed with the A form factors 
used in [l|. Dotted: RFG OB response. 



for the case of 12 C where the sizes of the correlation re- 
sponses, relative to the OB, are about 20% smaller than 
for 56 Fc. This is what one would expect, since the number 
of correlated pairs increases with A(A — l)/2. Moreover, 
since the estimated value of e depends on the nuclear 
radius, Eq. ([39]) indicates that one should use larger e- 
values for lighter nuclei, which in turn would reduce even 
more the size of the correlation responses. Thus we ex- 
pect an important A-dcpendcnce of correlations on the 
nuclear responses coming from the A-dcpendence of the 
escape width e. A more deailed study of this issue will 
be presented in forthcoming work. 



V. CONCLUSIONS 

In this work we have presented a fully relativistic model 
of inclusive two-particle emission reactions induced by 
electrons. Starting with the free relativistic Fermi gas 
we have considered all Fcynman diagrams in a perturba- 
tive expansion of the scattering amplitude with one pho- 
ton and one pion exchange producing 2p-2h excitations. 
Those diagrams can be classified in two sets, namely 
MEC and correlation currents. In the latter there is a 
nucleon propagator that can be put on shell giving a dou- 
ble pole from (jjq — E p + ie)~ 2 when taking the square 
of the current matrix element. The corresponding 2p-2h 
response function diverges as 1/e when e — > plus addi- 
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tional lne terms. Giving a physical meaning to e as the 
escape width of the nucleus, namely, twice the inverse 
of the nuclcon propagation time, the fact that the cor- 
responding response is infinite is related to the infinite 
extension of the Fermi gas. Using a finite value of e we 
account for the finite size of the nucleus, hence getting a 
finite result. Having no way to compute e in a Fermi gas, 
we take it as a parameter. Estimating in a crude way a 
value around ~ 200 MeV, we have made an exploratory 
study of the results as a function of e. The correlation ef- 
fects decrease with increasing e. Our analysis shows that 
the assumption e ~ 200 — 300 MeV is not unreasonable, 
whereas for smaller e- values the correlation contribution 
increases significantly in the QE region. 

Within this framework we have studied the properties 
and effects of the different 2p-2h contributions and other 
ingredients of the model on the transverse and longitu- 
dinal response functions of 56 Fe for intermediate to high 
momentum transfer. The MEC give rise to a wide peak 
in the region of the A resonance that dominates the T 
response. In the L channel the MEC are small for low 
momentum transfer, but they importantly increase for 
high momentum above the QE peak where their contribu- 
tion is of the same size as the OB longitudinal response. 
Concerning the correlations, they add to the MEC in the 
high-energy tail and are of the same order of magnitude. 
The contribution of the correlations is similar in the L 
and T responses. 

The main goal of this paper has been to study the ef- 
fect of 2p-2h pion correlations in the L and T response, 
analyzing the properties of these effects as a function of 
a single parameter e. In future work we plan to inves- 
tigate more physically founded ways to "fine tune" this 
parameter, including its dependence on kinematics and 
nuclear species. Finite-size calculations in conjunction 
with the use of semi-phcnomcnological fits of the nucleon 
spreading width or fits to existing (e, e') data will also be 
explored. 
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different channels PP (two protons), NN (two neutrons) 
and PN (proton-neutron) emission. 



1. PN channel 

We first consider the channel in which we eject a PN 
pair. In this case there is no symmetry in the wave func- 
tion and we assume that the first hole is a proton and the 
second is a neutron, i.e., the initial isospin wave function 
is | PA). The final state can be \PN) or \NP) depending 
on if there is or is not charge exchange. 

In the case of the MEC seagull and pion-in-flight, Fig. 
l(a-c), this is the only channel which contributes. The 
isospin operator is 



(Al) 



where repeated indices are meant to be summed. The 
relevant isospin matrix element is obtained by operating 
over a PN state 



(NP\U\PN) 



-2i. 



(A2) 



In the case of the correlation current we find four 
isospin operators for the diagrams of Fig. 2, including 
the isospin dependence in the single nucleon current T M , 
namely 



' a ' a J ' a ' a ' 

p,u(l) (1) (2) (1W2) (2) 



(A3) 
(A4) 



Operating over the initial \PN) state we obtain 



t^t^T^\PN) = 2T» N \NP)-T» N \PN) (A5) 

t WT^ 2) t^\PN) = 2T^ P \NP) -T^IPN) (A6) 

t^T^t^IPN) = 2T^ P \NP) -T^IPN) (A7) 

T^tPtW\PN) = 2T^ N \NP) -r" p |PA).(A8) 

In the case of the A current, diagrams of Fig. 1 (d-g), 
we find the following isospin operators 



T (l) T t(l) (2) r (l) T f(l) (2) 
(l) T (2) T t(2) T (l) T (2) T f(2) 



(A9) 
(A10) 



where Tj are the | — > ^ isospin transition operators ver- 
ifying 



(All) 



Appendix A: Isospin matrix elements 

In the model used in this work we compute explicitly 
the isospin matrix elements of the current operator in the 



For instance we have 

T (l) T t(l) T (2) ; 



T (l) T t(l) T (2) 
1 3 1 a T a 



H r (2) 

3 z 



T (D x T (2) 



T (D x T (2) 



(A12) 
(A13) 
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and similarly changing 1h2. Operating over the initial 
\PN) state we obtain 



T^T^r^\PN) = 



~\NP)-^\PN) (A14) 
1\NP)-^\PN) (A15) 



tPtWtP 2) \PN) = ~^\NP)+ Z -\PN) (A16) 
r (i) T f T t(2)| PiV) 



-\\NP) + ||P7V) . (A17) 



2. PP Channel 

In the case of two proton emission only the A and 
correlation diagrams contribute. In the case of the corre- 
lations, the isospin operators over the initial \PP( state 
give 



r ( a ] T^T^\PP) =P lP |PP) 



(A18) 



and exactly the same result for the remaining three op- 
erators. 

In the case of the A we have 

tPt1^tW\PP) = \\PP) (A19) 

and exactly the same answer for the remaining three op- 
erators. 



Appendix B: Integration of the energy delta function 

The 9-D integral for the 2p-2h response functions is of 
the type 

J d 3 p' 1 d 3 h 1 d 3 h 2 5{E 1 +E 2 +u J -E' 1 ~E 2 ) /(hl,h2,p' 1 ,p , 2 ) , 

(Bl) 

where p' 2 = hi + h 2 + q pi . The delta function allows us 
to perform one integration analytically imposing energy 
conservation. Therefore, for fixed values of the two hole 
momenta hi , h 2 and for fixed values of the two emission 
angles 8[, (/)[ of particle 1, we can integrate over the mo- 
mentum p' 1: fixing the energy of the first particle. To this 
end we change variables p[ —> E' = E[ + E' 2 . Taking into 
account that both energies E[ and E' 2 depend on p[ to 
compute the Jacobian of the transformation, we obtain 



dE' 



p'l 


Pa' P'l 


E'i 


E 'lV'l 



(B2) 



where the momentum of the final nuclcon for fixed emis- 
sion angles 6[ , 4>[ is obtained by solving the energy con- 
servation equation. This is a second degree equation with 
two solutions given explicitly by 



Pi 



V ± Vq\ 1 



(B3) 



3. NN Channel 

Once more only the A and correlation diagrams con- 
tribute. In the case of the correlations, we have 

t^t^T^INN) =T^ N \NN) (A20) 

and the same for the remaining three operators. 
Finally, for the A we have 

tWtZ {1) tM\NN) = ~\NN) (A21) 

and the same answer again for the remaining three oper- 
ators. 



where 





= V 2 


a 




b 


= E' 2 -p' 




= E' 


V 


= p' cos /3[ 



(B4) 

(B5) 
(B6) 
(B7) 



E' = Ei+E 2 +uj is the final total energy, p' = hi+h2+q 
is the final total momentum and P[ is the angle between 
p' 2 and p'. To compute the integral we add the contri- 
butions from these two solutions, corresponding to two 
possible final states compatible with energy-momentum 
conservation. 
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